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ABSTRACT 

Aims. SNR RX J 17 13.7-3946 is perhaps one of the best observed shell-type supernova remnants with emissions dominated by ener- 
getic particles accelerated near the shock front. The nature of the TeV emission, however, is an issue still open to investigation. 
Methods. We carried out a systematic study of four lepton models for the TeV emission with the Markov chain Monte Carlo method. 
Results. It is shown that current data already give good constraints on the model parameters. Two commonly used parametric models 
do not appear to fit the observed radio. X-ray, and y-ray spectra. Models motivated by diffusive shock acceleration and by stochastic 
acceleration by compressive waves in the shock downstream give comparably good fits. The former has a sharper spectral cutoff in 
the hard X-ray band than the latter Future observations with the HXMT and NuSTAR may distinguish these two models. 
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1. Introduction 

Energetic particles produce radiations over a broad energy range 
under typical astrophysical circumstances. To study the ener- 
getic particle population in high-energy astrophysical sources, 
observations over a broad spectral range are needed. Multi- 
wavelength observations of individual sources usually give 
sparse spectral data points, which can be fitted reasonably well 
with a simple power-law or broken power-law function. A fea- 
tureless power-law distribution implies the lack of characteristic 
scales or distinct processes in the system. The nature of particle 
acceleration processes are still a matter of debate nearly a cen- 
tury a fter the disc overy of high-energy particles from the outer 
space (lButtll2009l) . 

Recent progress in observations of high-energy astrophysi- 
cal sources has resulted in more detailed spectral coverage, and 
spectral features start to emerge. These features are produced 
by the underlying physical processes with well-defined charac- 
teristic scales in space, time, and/or energy, and they play es- 
sential roles in advancing our under standing of the acceleration 
mechanism (Berezhko & Krvmskvlfl988) . Diffusive astrophys- 
ical shocks are considered as one of the most important parti - 
cle accelerators (!Topt vgin|[T980l;lDrurvlll983l:lBerezhkolll996h . 
The acceleration of particles by compressive motions of astro - 
physical flows is also a very generic process (|Ptuskinlll988h . 
Both processes may be important in the acceleration of particles 
dominating the radiative cha racteristics of the shock associated 
with SNR RX J1713.7-394 (iFan. Liu & FrverllMO ). Detailed 
X-ray and y-ray observations of this source have revealed several 
features: 1) the X-ray and y-ray emissions ar e well correlated 
in space (lAharonian et al. 1120061 lPlagall2008h : 2) the emission 
spectrum sh ows a clear high-energy cutoff in both the X-ray and 
TeV bands (" Aharonian et al. 2007), and the cutoff in hard X- 
rays appears to be sharp (Tanaka et al. 2008), implying a sharp 
cutoff in the electron distribution producing the observed radio 
to X-ray spectrum through the synchrotron process; 3) the y-ray 



spectrum has a broad convex shape ( 'Funkll2b09h : 4) there are 
bright X-ray filaments with a width of ~ 0.1 hght y e ar varying 
on a timescale of about a year dUchivama et al. l2007h . lLiu et al. I 
(2008) showed that these emissions may be attributed to a single 
population of electrons with the y-rays produced via the inverse 
Compton scattering of the background photons by relativistic 
electrons. 

In light of recent high-resolution observations with Suzaku 
and Fermi, we use the Markov chain Monte Carlo (MCMC) 
method to constrain parameters of four possible lepton models. 
It is shown that simple models with a power-law electron dis- 
tribution, which cuts off at TeV energies, give a poor fit to ei- 
ther the X-ray or y-ray spectrum. Motivated by the mechanism 
of diffusive shock acceleration (DSA) and stochastic accelera- 
tion (SA) by compressive motions, we introduce two alterna- 
tive models, both of which give acceptable fits to the observed 
spectra (Section|2]l. In Section|3] we discuss the implications of 
these results and future work needed to improve these models. 
Conclusions are drawn in Section |4] 

2. Spectral fits with lepton models 

In the lepton scenario, the radio to X-ray emissions are produced 
through the synchrotron process of relativistic electrons, and the 
y-ray emission is produced through the inverse Compton scatter- 
ing of the background radiation fields by the same electron popu- 
lation. If the background magnetic fields are relatively uniform, 
similar to the background radiation fields, the model can nat- 
urally explain the spatial correlation between the X-rays and y- 
rays (|Plagall2008l) . As pointed out bv lLiu et al. 1 (2008). the mag- 
netic field required to produce the X-ray flux agrees with what 
is required to produce an X-ray spectral cuto fl^ by an electron 
population cutting off in the TeV energy range (lAharonian et al. I 
2007; Ellision et al. 2010). 

There are at least four parameters in a specific model. In 
this work we employ an MCMC technique suitable for multi- 
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parameter determination to search all the model parameters. The 
Metropolis-Hastings sampling algorithm is adopted when deter- 
mining the ju mp probabiHty fro m one point to the next in the pa- 
rameter space (lMacKavll2003h . The MCMC approach, which is 
based on the Bayesian statistics, is superior to the grid approach 
with a more efficient sampling of the parameter space of inter- 
est, especially for high dimensions. The algorithm ensures that 
the probability density functions (PDF) of model parameters can 
be asymptotically approached with the number density of sam- 
ples. Starting with an initial parameter set Pq, which should lie 
in a reasonable parameter space, one obtains the likelihood func- 
tion of Pq: -Co(Po) '^ exp(-;t'^(Po)/2|. Then another parameter 
set Pi is randomly selected with the corresponding likelihood 
function £,\. This parameter set is accepted with the probability 
a - min{ 1 , Xi /-Co}- If the new point is accepted, it becomes the 
starting point for the next step, otherwise, we reset Pi = Po. This 
procedure is repeated to derive a parameter chain that determines 
the PDF of model parameters. More details of this procedure can 
be found in Neil (1993), Gamerman (1997), andMacKay (2003). 
I: Without detailed modeling of the particle acceleration pro- 
cesses, a power-law electron distribution with an exponential 
high-energy cutoff is often assumed to fit the observed spectrum: 



N{y) ocy P exp(-y/yc) 



(1) 



where p is the spectral index, y^ the Lorentz factor of the high- 
energy cutoff Ec - jcineC^, nie and c are the electron mass and 
the speed of light, respectively. The integration of N{y) over the 
Lorentz factor y gives the particle number density. It has been 
shown that this model gives a sufficient fit to previous obser- 
vations of SNR RX J1 7 13.7-3946 when t he interstellar radia- 
tion field is considered (iPorter et al. Il200a) : however, more de- 
tailed studies suggest that such a model systematically underes- 
timates the photon flux at the low-energy end of HESS obser - 
vations (iTanaka et al. I l2008t IZirakashvili & Aharoniaii] 1201 Oh . 
In this paper, \ye adop t the interstellar radiation field spectrum 
of Porte ret al.l ([2006"). There are then four model parameters. 
Besides the magnetic field B, p, and jc, we need a normalization 
for the relativistic electrons /Stot, which is the total energy of elec- 
trons with y > 10. Figures[T]and|2]give the best fit of this model 
to the data and the probability density of the model parameters. 
All four model parameters are well constrained. However, de- 
tailed examination of residuals of the spectral fit shows that there 
are several data points below 1 and near 10 TeV with normalized 
residuals significantly greater than 3. 

n: |Liu et al. 1 12008) suggested that the fit to the y-ray spec- 
trum can be improved by considering a more gradual cutoff in 
the electron distribution: 



N{y) ocy P exp -{ylycf^^ . 



(2) 



Several mechanisms can cause 

of the high-energy cutoff (iBerez hko & Krvmskv 



Becker. Le & Dermerl 120061 : IZirakashvili & Aharonian 



variations in the shape 
1 11988 : 
l] 



2007; 



Blasi 2010). However, giv en the sharp cutoff of the X-ray spec- 
trum, iTanaka et al. I (l2008l) claim that the electron distribution 
must have a very steep cutoff. Indeed, by using their data and 
performing the MCMC fit, we obtain radio fluxes nearly one 
order of magnitude below the observed values. The model also 
significantly overestimates the hard X-ray fluxes. To increase 
the weight of the radio and y-ray data in the model fit, we 
artificially increase the errors of the X-ray data from Suzaku 



by a factor of 2 Q (Figs. \T\ and |2]l. The normalized residuals 
are indeed improved significantly with a reduced x^ of 1.2, 
in agreement with Liu et al. (2008), where cruder X-ray data 
fro m ASCA obse rvations were used. However, as emphasized 
bv .Tanaka et al. I (12008,) . the sharp cutoff in the X-ray spectrum 
is significant and th e error estimate of the Suzaku observation 
appears to be robust (Takahashi e t al. 1120081) . This simple model 
therefore cannot fit the X-ray spectrum. 

Ill: The particle distribution is determined by the injec- 
tion process and the spatial diffusion coefficient x in the DSA 
model. The time-dependent electron dis tribution with a constant 
injection rate may be approximated as (Fo rman & Drurv 11 19831 : 
Berezhko & Krvmskv ..1988;.Drurv ..1991.) 



N(y) <x y~''exp[-9;r(r)/t/'Tiife] , 



(3) 



where Tufe is the supernova lifetime since the onset of particle 
acceleration and U the shock speed. A very sharp high-energy 
cutoff can be produced when gyro-radii of high-energy particles 
exceed the characteristic length of the magnetic field in the back- 
ground plasma Ij. Particles at such high energies do not perform 
gyro-motions around the chaotic magnetic field and instead in- 
teract randomly with the background fields. The spatial diffusion 
coefficient of these particles increases as the square of the rela- 
tivistic Lorentz factor in the absence of waves on scales greater 
than Id- Consequently, they can readily cross the acceleration site 
without interactions. We assume a spatial diffusion coefficient 



^^^' ' ' y2m,cV(3eByo) 



for 
for 



y <yc 
y>yc 



(4) 



where e is the elementary charge, 77 > 1 depends on the turbu- 
lence intensity, fi - \j2 and 1/3 for particle interactions with a 
spectrum of p l asma waves following the Kraichnan ( 1965) and 
iKoknogorov I (Il94lh phenomenology for the turbulence spec- 
tral evolution, respectively. The characteristic length of the mag- 
netic field is then given by the gyro-radius of electrons at y^: 
Id - ycineC^leB. 

Compared with the previous two models, an extra pa- 
rameter is introduced here: ri/Tuf^U^, which is proportional 
to the electron acceleration timescale. Figures [1] and |2] show 
the results with /3 = 1/2. The probability density of 
77(1.6kyr/riife)(4Mm s"'/t^)^ peaks near a typical value of 12. 
Since 77 has to be greater than 1, we have Tm^U^ > 0.02c 
lyr The gyro-radius of electrons at the cutoff' Lorentz factor 
yc ~ 2 X 10^ is ~0.01 pc comparable to the width of the ob- 
served highly variable X-ray filaments. The X-ray variability 
may then be attributed to spatial diffusion of X-ray emitting elec- 
trons from a high-density region formed presumably via an inter- 
mittent process. It is interesting to note that the best fit spectral 
index p - 1 .92 is less than 2, implying acceleration by multiple 
shock s or strong nonlinear effects ( Zirakashvili & Aharonian] 
l2010h . The model with/? =1/3 gives a similar fit with a slightly 
higher value of 77 and sharper har d X-ray cutoff. 

IV: iFan. Liu. & Frverl (l2010l) have carried out detailed mod- 
eling of electron acceleration by compressive motions in the 
shock downstream and argued that it might also explain these 
observations. Here we follow the same treatment and also con- 
sider the effect of incompressive motions, which can dominate 
the spatial diffusion and therefore the escape of energetic elec- 



' Given the complexity of the instrumental background and calibra- 
tion, and unknowns, such as the estimated cosmic X-ray background, 
such an artificial increase may well be reasonable. 
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tro ns from the acceleration re gion. The acceleration rate is given 
bv (iBvkov & Toptvgin|[T99l 



Dyyly^ 



4nx r'' 
9 Jl„ 



dk- 



k'^S(k) 



(5) 



where x, K vf, and S are the spatial diffusion coefficient, the 
wave-number, speed, and intensity of compressive waves, re- 
spectively. The overall wave intensity J S {k)'Xnk^dk — ^vj, with 
^ < 1 for subsonic turbulence. Q The turbulence is generated 
on a scale L near the shock front, and Id corresponds to the 
scale where the turbulence speed is comparable to the Alfven 
speed va- Relativistic electrons with a gyro-radius less than Id 
have a mean-free path comparable to Id since waves on smaller 
scales are subject to strong transit time damping by thermal 
background ions. We then have the diffusion coefficient 






for y < ydx) , 
for y > ydx) .. 



(6) 



where x indicates the spatial location in the downstream, /,/ - 
ycnigC^ jeB, and we have assumed Bohm diffusion for y > yc- 
Due to gradual decay of the turbulence intensity, related quanti- 
ties evolve with x in the downstream. Equation (|6]l is quite dif- 
ferent from equation Q, where particle scattering by the back- 
ground turbulence is parameterized without considering details 
of the turbulence structure. Considering the presence of strong 
turbulence with m > va between the scales of U and L, Bohm 
diffusion may be appropriate, and as we show below, faster dif- 
fusion may not be able to accelerate electrons efficiently to ac- 
count for the observations. The intensity of the incompressive 
motions T{k) can be determined from the energy conservation: 
J[2r(^)H-5 (K)\\Tik^dk — 3u^, where 3m^/2 is the total turbulence 
energy density per unit mass. Further downstream, where Vf is 
high (> 6'^^m), the incompressive motions vanish completely. 
To derive the escape time T^sc of electrons from the acceleration 
region, we need to consider the enhancement of the spatial dif- 
fusion c oefficient by ffie turbulence motions. According to Eq. 
(4.30) in iBvkov & Toptygin I (Il993h . ffie effective diffusion co- 
efficient ;t', can be derived from 



X* 



8n r 
X+^^ 
^X* Jk 

3 Jk„ 



kd 



T(k)dk + 

1 2(kY,-vl) 



(7) 



[k^xUvj (k^xl + vD^l 



k^S(k)dk, 



where we have ignored the decay of incompressive motions. 
Then we have 



Tesc = L Ix* 



(8) 



With the above treatment, we can study particle acceleration 
by compressive motions all the way to the shock front , whic h 
is different from the case studied by ^ Fan. Liu. & Fryer] (l2010h . 
where only the subsonic phase with v/r > m is considered. Near 
the shock front, we assume ^ - 1 /2 so that equation (|5]l is ap- 
phcable. Since ffie comp ressive wave intensity is lower than in 
iFan. Liu. & Fryer] (1201 Oh and the acceleration timescale of TeV 
electrons becomes comparable to the supernova lifetime, a time- 
dependent treatment is necessary. Boffi ffie acceleration and 



^ Here we adopt the formula given by iBvkov & Toptygin lll993l). 
which gives an acceleration rate a factor of 2 lower than that of lPtuskin] 
CHS). 

^ At the point, where vp = u, the compressive wa ve intensity used 
here is 3 times lower than in lFan. Liu. & Frverl ( l2O10h . 




1C 


10 


10 


10 10 10 10 


10 




10 


10 10 


1C 


'' 


2 


2.3 


--; 


-^^-■Sf-W^¥^^^'^-'?---~-'.7. 


— T 


-^ 


-_-_ 


— ^ ^ .: 


='-\ 


°2 


b 


1.2 


--; 


-^^sf**^***'*""^'-^^;^; 


r--_ 


- 




1 


- 


-2 










b 
-2 


1.7 


--; 


\'^^Hf^!n^f^^'f^-'-A^ 


■* — - 


:--^ 


^--' 


"*— -'-:-•"--■-:--; 


='--_ 


-2 


2 
4 


1.9 


--; 


;^^^#^?s^--r^s 


— -_ 


-^ 


_-_- 


"^— -'-:'^v.-,_; 


--: 






10 10 10 10 10 10 



Energy ( eV ) 



Fig. 1. Model ffi to the spectrum o f SNR RX J1713.7-3946. The 
TeV data is fr om PUiaronian et al. I (l2007i) . t he Fermi data i s from 
iFunkI (l2009h . and the X-ray data is from iTanaka et al. I (l2008h 
(Thanks to Jun Fang). The thick solid, dashed line and the thin 
solid, dashed line in the upper panel correspond to the diffu- 
sive shock, stochastic acceleration model, and the model with 
an exponential, more gradual cutoff in the electron distribution, 
respectively. The low and high-energy spectral hump are pro- 
duced by relativistic electrons through the synchrotron and in- 
verse Comptonization process, respectively. The lower panels 
show the normalized residuals. From top to bottom, they are for 
the model with an exponential and more gradual cutoff, the DS A, 
and SA model. For the model with a gradual cutoff, the X-ray er- 
rors have been artificially increased by a factor of 2 to improve 
the fit. The reduced x^ of these fits are indicated in the top-left 
corner See the text for details. The instrumental sensitivity of 
NuSTAR and HXMT are indicated by the dotted and dot-dashed 
lines, respectively. 



escape timescale explicitly depend on the diffusion coefficient 
X, which is independent of y for y < y^. Assuming that elec- 
trons are injected at the Lorentz factor y,„, = 10 with a constant 
rate, ffie time-dependent results o f Becker. Le & Dermer I (l2006t 
see also lCowsik & Sarkarll984l:rPark & Petrosionlll995h can be 
used to describe the evolution of the electron distribution func- 
tion at energies below y^. Near the cutoff energy y^-, the nontriv- 
ial dependence of the acceleration and escape timescales on the 
Lorentz factor demand numerical solutions. We assume the time- 
dependent distribution function has a high-energy cutoff given 
by exp[-(Tfl(./7"cic)'''^] (For details, see Fan, Liu, & Fryer 20l3)- 
Figures ^ and (|2| show the results for U - 4000 km/s with the 
Kraichnan phenomenology. The quality of ffiis fit is comparable 
to the DSA model. T he deduced model pararn eters also agree 
with those obtained in lFan. Liu. & Fryer] (|2010|) . 



3. Discussion 

The SA model (Fig. [U predicts a higher hard X-ray flux than 
the DSA model since the diffusion coefficient of the DSA model 
has a stronger dependence on the electron energies leading to 
a sharper high-energy cutoff. The electron distribution in the 
SA model includes contributions from broad regions in the 
downstream, which also makes the overall electron distribution 
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L (I.OelScm) 



Fig. 2. Probability density function of model parameters normal- 
ized at the peak value. Different line types correspond to differ- 
ent models as explained in Figure [T] and the text. 



broader Future observations with NuSTAR and HXMT may be 
able to distinguish these two models. 

In actual turbulent astrophysical shocks, both the DSA 
and S A mechanisms may contribute to the electro n acceler- 
ation (AchterbergJ 119901; iBvkov & Toptv^in 1 11993). The dis- 
tinction between the two res ides in the dissipation structure 
(iBerezhko & Krvmskvlll988h . Supersonic shock fronts domi- 
nate the acceleration in the DSA model. The SA model domi- 
nates in the subsonic phase. It is possible that the shock down- 
stream of SNR RX 1713.7-3946 has both a supersonic and sub- 
sonic phase turbulence with the former closer to the shock front. 
The electron acceleration is therefore a continuous process in 
the turbulent downstream. In the S A model of this paper, we in- 
tentionally remove the DSA by requiring that the compressive 
waves be subsonic in the downstream. A self-consistent treat- 
ment of the turbulence evolution in the shock downstream is 
needed for particle acceleration near astrophysical shock fronts. 

All the models appear to systematically underproduce y- 
ray flux in the Fermi energy band and overproduce flux near 
10 TeV. The latter could be caused by inhomogeneities of the 
background magnetic field, which implies lower cutoff energy 
of the electro n distribution and less high-energy y-ray flux 
dTanaka et al. 1 12008 ). The former requires a broader electron 
distribution from the GeV to TeV energy range. A combina- 
tion of DSA and SA model with the former dominating the 
higher energy particle acceleration and the latter enhancing the 
lower energy ones may address this issue. Such a scenario is 
possible if both supersonic and subsonic turbulence are pro- 
duced by the supernova shock (IBerezhko & Krvmskv I Il988t 
iBvkov & Toptvgin1ll993h . Contributions to GeV-TeV emission 
from decay of neutral pions produced by inelastic collisions of 
relativistic protons with the background ions may also e xplain 
thes e residuals. How e ver, a s shown by lEUision et al. I (l2010l) 
and lKatz & Waxmanl (l2008h . relativ istic protons can not be the 
dominant TeV emission component dZirakashvili & Aharonian I 
12010 "). Including this component is beyond the scope of the cur- 
rent investigation but may give a good constraint on the rela- 
tive acceleration efficien cy of relativistic electrons and protons 
(lKatz&Waxmanll2008h . 



4. Conclusions 

With recent high spectral-resolution observations of SNR RX 
1713.7-3946 in X-rays and y-rays, we show that, in general, 
models with simple parametric descriptions of energetic particle 
distribution give a poor fit to the broadband spectrum. Details of 
relativistic electron acceleration may be probed with advanced 
models. Although both the DSA and SA model give acceptable 
fits to the spectrum, the SA model predicts a higher hard X-ray 
flux than the DSA model, which may be tested by future ob- 
servations. The spatial distribution of energetic electrons in the 
SA and DSA models is also different, and the source structure 
can be used to distinguish these two models as well. The results 
agree with the scenario where energetic electrons are accelerated 
by interacting diffusively with a turbulent electromagnetic field 
and produce nearly all of the observed emissions from the shell 
of SNR RX 1713.7-3946. Turbulence evolution at astrophysical 
shocks appears to play a key role in advancing our understanding 
of particle acceleration processes. 
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